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Abstract 



We develop a simple systematic method, valid for all strengths of disorder, 
to obtain analytically for the first time the full distribution of conductance 
P(g) for a quasi one dimensional wire in the absence of electron-electron 
interaction. We show that in the crossover region between the metallic and 
insulating regimes, P{g) is highly asymmetric, given by an essentially 'one 
sided' log-normal distribution. For larger disorder, the tail of the log-normal 
distribution for g > 1 is cut off by a Gaussian. 
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Since the discovery of the absence of self averaging in mesoscopic disordered systems 
0, the study of the full distribution of conductance has attracted a lot of attention [Q. In 
particular, while the metallic regime is well described by a Gaussian distribution, the mo- 
ments of the conductance fluctuations become of the same order of magnitude as the average 
conductance on approaching the localized regime. In such cases the average value becomes 
insufficient in describing properties of disordered conductors and the full distribution must 
be considered. Recently, numerical support for the existence of a new universal distribution 
at the metal-insulator transition a broad distribution of the critical conductance at the 
integer quantum Hall transition || , as well as the expected multifractal properties associated 
with the critical regime |J have increased the interest in the conductance distribution in 
the intermediate regime, between the well studied universal conductance fluctuations in the 
metallic limit and the log-normal distribution in the deeply insulating limit. However, even 
for a quasi one dimensional (Id) system where there is only a smooth crossover between 
the metallic and insulating regimes, there is no analytic result available for the conduc- 
tance distribution in the crossover regime. So far only the first two moments have been 
obtained for all strengths of disorder 0, using the Id supersymmetric nonlinear a model 
0. This model has been shown || to be equivalent, in the thick wire or quasi Id limit, to 
the Dorokov-Mello-Pereyra-Kumar (DMPK) equation [[J which describes the evolution of 
the distribution of the transmission eigenvalues with increasing wire length. 

In this work we develop a simple systematic method to evaluate directly the full distri- 
bution of conductance for a thick quasi Id wire (mean free path I 3> the width), starting 
from the solution of the DMPK equation. The main result of the paper is that although 
there is no phase transition in quasi one dimension, the crossover region between metallic 
and insulating regimes is highly non trivial, and shows a remarkable 'one-sided' log-normal 
distribution. Recent numerical studies of a quasi Id system in the quantum Hall regime 
have shown highly asymmetric log-normal distributions in the crossover region . We ex- 
pect similar qualitative features to exist in the critical regimes in higher dimensions as well. 
Indeed, numerical studies near the integer quantum Hall transition in two dimensions as 
well as the Anderson transition in three dimensions also point to asymmetric distributions 
of the critical conductance jnj. In addition, we predict that even the insulating regime 
should have a sharp cutoff in its log-normal tail near the (dimensionless) conductance g ~ 1. 
In particular, we show that the conductance distribution in the insulating regime (in the 
absence of time reversal symmetry) has the form 

{/ xi sinh 2xi „— ra? „ ^ i . 
V i-* e , ' 9<L ' (1) 
yfige-*- 1 ?, g>l 

where x\ = cosh" 1 (l/ v /^) and the parameter V = £/L, where £ = Nl is the quasi Id 
localization length, N is the number of transmission channels, and L ^> I is the length of 
the conductor. The parameter a is the value of F" given in (10) evaluated at X2 = 2/7rr, 
and tends to |exp[8/7rr] for T <C 1 in insulators. Note that for g <C 1, £i ~ \ia(2/^/g), 
and the distribution is log-normal, centered at — hig = 1/T. However, for g > 1, the tail is 
cut off by an exponential function over an exponentially narrow scale in 1/T. The results 
for two different values of T, 0.7 and 0.2, are plotted in fig. 1. The main point is that for 
very strong disorder, the typical value of g is much smaller than unity, so the peak of the 
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distribution P[ln(g)] is very far away from g <~ 1. In this case the exponential cutoff at 
g ro 1 is less relevant. However, even for large disorder, a sharp cutoff in the tail for g > 1 
always exists, as shown for T = 0.2 in fig. 1. At intermediate strength of disorder, the peak 
of the distribution is close to the cutoff, and the distribution becomes highly asymmetric. In 
particular near the crossover between metallic and insulating behavior, the peak is at g ~ 1, 
and we obtain a 'one-sided' log-normal distribution, as shown for V = 0.7 in figure 1. 

As a check of the scope and validity of the method developed here, we obtain the exact 
universal conductance fluctuation in the weakly disordered metallic regime with the expected 
Gaussian distribution as well as the correct mean and variance of the log-normal distribution 
in the strongly disordered localized regime within the same unified framework. Systematic 
corrections as a function of disorder can be obtained from both the metallic and the insulating 
limits. Note however that the analytic results presented here near the crossover regime are 
only semi-quantitative, due to the approximate analytical evaluations of certain integrals. 
A more quantitative result is possible based on numerical evaluations of these integrals. 

We first briefly outline the method. The probability distribution p(X) of the N variables 
Aj, related to the transmission eigenvalues Tj of an N-channel quasi Id wire by A» = (1 — 
Ti)/Ti, satisfy the well known DMPK equation ||. The solution of this equation can be 
written in the general form |TI| p(X) oc exp[— j3H(\)], where H(X) may be interpreted as 
the Hamiltonian function of N classical charges at positions A«. The symmetry parameter 
(3=1,2 or 4 depending on the symmetry of the ensemble 0. The Hamiltonian depends on 
the parameters L, N and / only in the combination T = Nl/L. We will consider the quasi Id 
limit where both N and L approach infinity keeping T fixed. The dimensionless conductance 
in terms of Aj is given by g 
be written as 



J2i y+a~ EH- ^he distributon of conductance can therefore 



P(g) 



dr 

Z J-oo 2% Jo 
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+ A, 



(2) 



where Z is a normalization factor. In the metallic regime g ^> 1, the Aj are all very close 
to each other so that a continuum description can be used with a density of A finite be- 
tween zero and an upper cutoff given by the normalization condition. This approximation 
describes the universal conductance fluctuations in the metallic regime fl4 |. In the deeply 
insulating regime on the other hand, all A» are exponentially large and separated exponen- 
tially from each other, and the conductance is dominated by the smallest eigenvalue. This 
approximation describes the log-normal distribution in the deeply insulating regime JT5|. It 
is clear however that none of the above descriptions can be used in the crossover regime, 
where the smallest eigenvalue is neither zero, nor exponentially large. Nevertheless, it turns 
out that it is possible to combine the essential features of the two descriptions and develop 
a simple and systematic procedure to study the conductance distribution at intermediate 
regimes. For simplicity, we will discuss the case [3 = 2 only. 
The basic idea is the following: 

1) We first separate out the lowest eigenvalue Ai and treat the rest as a continuum with 
a lower bound at A2 > Ai. Note that this approximation can be systematically improved by 
separating out the lowest n > 1 eigenvalues and treating the rest as a continuum. 

2) The continuum part can be written as a functional integration on the generalized 
density p(A), and the distribution (2) can be rewritten as 
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P{g) = z ^ lT9 I dX ± J dX * J D ^ exp[-F(A!, A 2 ; p(A); r)]. (3) 



Here the 'Free energy' 

P(A) 



F(X 1 , A 2 ; p(A); r) = PH(X 1 , A 2 ; p(A)) + *r 



1 + Ai A 2 1 



+ A 



(4) 



contains the 'edge' separating out Ai as well as the 'source' terms proportional to r, plus 
the continuum version of the Hamiltonian of the form H = £^ w(A i; Xj) + £f ^(A;). The 
upper limit 6 is given by the number conservation dXp(X) = N — 1. 

3) We obtain the density by minimizing the Free energy with respect to p(A), keeping 
Ai and A 2 fixed. This gives 



roo 

Jo 



dX'u(X + A 2 , A' + A 2 )p sp (A + A 2 ) = 2V tot (X + A 2 ), (5) 

where we have shifted the lower limit to zero, and V tot (X) = V(X) + + u(Xi, A). After 
taking a derivative on both sides of (5), the kernel can be inverted to obtain the saddle point 
density p sp (X). 

4) From the density, we obtain the saddle point Free energy 
P P f b i \ t t" ,n ,n , m ,M ^ . ir 



sp 



2 dXV tot (X) Psp (X) + (3V(X 1 ) + (6) 



5) Since V to t and therefore p sp are both linear in r, the Free energy is quadratic in r and 
can be written in the form F sp = F° + (ir)F' + ^-F". The integral over r in eq. (3) can 
then be done exactly. The result is 

1 roo roo ( n 

6) At this point, the integrals over Ai and A 2 can be evaluated numerically. Instead, we 
use saddle point approximation to do the integrals in order to obtain an analytic expression 
for P(g). Solving for = for i — 1,2 to determine the saddlepoint values of Ai and A 2 , 
we obtain the distribution as a function of the conductance g, in terms of the parameter T. 

In the above approach, if we set both Ai and A 2 equal to zero, we obtain the correct 
universal value 2/15/3 for the variance of g. This is consistent with the picture that in 
the metallic regime, the eigenvalue density can be treated as a continuum from zero to an 
upper cut off b. As disorder is increased beyond the metallic regime, this picture starts to 
break down. In particular, the smallest eigenvalue is pushed to a finite distance from zero 
depending on the strength T, so that the continuum picture at the edge no longer holds. 
The correction to the metallic behavior is captured in the present approach by evaluating 
the shifts in Ai and A 2 within a variational scheme. Since the insulating regime is dominated 
by the smallest eigenvalue, this approach clearly captures the correct insulating behavior. In 
the crossover regime, both the separation of the smallest eigenvalue as well as the rest of the 
continuum become important. Note that if more accuracy is needed, one can in principle 
separate out more than one eigenvalue. 
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We now give some details. From the exact solution of DMPK eqn., the two and one 
body terms in the Hamiltonian of (2) are known to be [|l(J 

u(X, A') = -~ In | (A - A')(* 2 (A) - * 2 (A'))|; V(X) = ^r 2 (A), (8) 

where x(A) = sinfT 1 y/X. Note that the difference Au = u(X + A 2 , A' + A 2 ) — u(X, A') 
is negligible in the insulating regime and is a small correction proportional to A 2 in the 
metallic regime. Therefore to a first approximation, the shifted kernel u(X + A 2 , A' + A 2 ) 
can be replaced by the unshifted kernel w(A, A')- One can then use the saddle point density 
obtained from the unshifted kernel to calculate the correction due to the change in the 
kernel from the shift. This can be rewritten as an additional term A 2 V 2 + V to t = V e ff with 
an unshifted kernel in eq. (5). The unshifted kernel can then be inverted to give the saddle 
point density 

1 f°° d 
Psp (\ + A 2 ) = j ^ d\'K-\x{\) - x (X'))-^V eff {\X'\ + A 2 ), (9) 

where the inverse of the unshifted kernel is K~ l (t) = — (1/27T 2 ) / °° dqsin(qt)(l — e _?r9 ). The 
condition p sp (A) > for all A in eq. (9) requires A 2 — Ai > A c = (2/T7r) 2 . The free energy 
can now be obtained from eq. (6). The integrals for F" can be done exactly. In terms of 
the variables x\, x 2 , defined as sinh 2 ^ = Xi,i = 1,2, we get 



1 

~Z + 



3 ' 4x| sinh 2 2x 2 _ 

The integrals for F' and F° can be evaluated analytically in two limits. For X2 C 1, 

F'^T-b 1 Txl + -Jxl-xl 



(10) 



7T 3 



3tt 2 



F° ss — r 2 ^ - 2TxY\jxl + - ln(x 2 - x\) - h\x x , (11) 

where b\ ~ 0.89. In the other limit x 2 ^> 1, 

1 11 

F' P3 g — ; F ° ~ Tx i ~ - M x i sinh2xi) + -r 2 ^ - Tx 2 2 + x 2 . (12) 

cosh xi 2 3 

In the metallic regime, T ^> 1, and x 2 can be very small. Then eqs (10) and (11) give the 
correct mean conductance (g) = T and variance var(g) = 1/15. When T < 1, A 2 — Ai > A c 
requires a; 2 ^> 1. In this case the limit X\ ^> 1 corresponds to the insulating limit, but 
the limit x\ 1 corresponds to the intermediate case close to the crossover regime. We 
therefore study this regime within a saddle point approximation for the integrals (7). 

The condition J^- = has the solution coshxf = -4=, while the condition = has 

axi 1 y/g' dx 2 

the solution x 2 p = 1/r. This leads to the saddle point result S sp , to which the contributions 
from the fluctuations S^ 1 = In |<9 2 Sy<9x 2 | have to be added, leading to eq. (1). In the 
deeply insulating regime r<l, the above expression leads to the known mean and variance 
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(ln(l/g)) = var\\.n(g)]/2 = 1/T. However, note that since coshxi > 1, the saddle point 
solution exists only for g < 1. For g > 1, the X\ and £2 integrals are dominated by the 
boundary values at X\ = and X2 = 2/71T, which has been incorporated in (1). In fig. 2 
we show P(g) as obtained from eq. (1) for T = 0.7. The skewed shape and the exponential 
drop at g ~ 1 compare well with numerical results of |L1| for a slightly smaller value of 
T = 0.5. The difference in the T values simulates to some extent the correction terms to 
eq. (1) expected for values of T approaching unity. Also shown in fig. 2 is the result of 
a numerical integration of eq. (7) using eq. (11) for T = 1.6 in the metallic regime. The 
Gaussian shape of P{g) obtained for this rather small value of T is in very good agreement 
with the results of [ p"T|j . 

We note that according to the relation g = ji^i an y non-negligible P(g ~ 1) comes 
from the possibility that the smallest eigenvalue Ai can be close to the origin. However, given 
Ai <C 1, the logarithmic repulsion between eigenvalues generated from (8) forces the rest of 
the eigenvalues exponentially far when r < 1, so P(g > 1) is cutoff sharply. Since these 
arguments are quite general, we expect qualitatively similar features in higher dimensions as 
well, which should have important consequences for the universal conductance distribution 
in the critical regime. 

To summarize, we calculated the distribution of conductances P(g) for a quasi Id dis- 
ordered system using known results for the DMPK equation. In this case P{g) depends 
only on one parameter V = £/L, where £ is the localization length. In the crossover regime 
£/L ~ 1, we find that P(g) is given by a 'one-sided' log-normal distribution, cut off by a 
Gaussian tail on the metallic side (g > 1). We believe that this behavior could be generic 
for P(g) in the transition regime even in higher dimensions, provided the average of g at the 
transition or crossover region is of order unity. Our results can not be directly compared 
to the work of |]TB[ in d = 2 + e dimensions, because for e 1, (g) = 1/e is large, and 
the bulk of P(g) is located deep in the metallic regime. As proposed in | 17[| , the center 
of P(g) is then Gaussian, with power law tails oc g~ 2 ^ e . The latter results are peculiar to 
the behavior in 2 + e dimensions. One should keep in mind that the DMPK approach does 
not contain the effects of wave function correlations in the transverse direction, which are 
expected to be important in higher dimensions. Nonetheless, the similarities of the shape 
of P(g) in the crossover regime obtained here with the numerically determined P{g) in 3d 
at the critical point |12| appears to suggest that the generic behavior of P{g) is that of a 
log- normal distribution for g < 1 combined with a Gaussian cut off for g > 1. 

We are grateful to A. Mirlin for stimulating discussions and bringing refs. |TT| and [O to 
our attention. We also thank M. Fogelstroem for his help regarding numerical evaluations. 
This work has been supported in part by SFB 195 der Deutschen Forschungsgemeinschaft. 



6 



REFERENCES 



B.L. Altshuler, JETP Lett. 41, 648 (1985); P. A. Lee and A.D. Stone, Phys. Rev. Lett. 
55, 1622 (1985). 

For reviews, see Mesoscopic Phenomena in Solids ed. B.L. Altshuler, P.A. Lee and R.A. 
Webb, Elsevier, Amsterdam 1991; C.W.J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

B. I. Shklovskii et al, Phys. Rev. B47, 11487 (1993); P. Markos and B. Kramer, Phil. 
Mag. B68, 357 (1993). 

Z. Wang, B. Jovanovic and D-H. Lee, Phys. Rev. Lett. 77, 4426 (1996); S. Cho and 
M.P.A. Fisher, Phys. Rev. B55, 1637 (1997); D.H. Cobden and E. Kogan, Phys. Rev. 
B54, R17316 (1997). 

For a review, see M. Janssen, Phys. Rep. 295, 1 (1998). 

M. R. Zirnbauer, Phys. Rev. Lett. 69, 1584 (1992); A. D. Mirlin, A. Muller-Groeling 
and M. Zirnbauer, Ann. Phys. (N.Y) 236, 325 (1994). 
K. B. Efetov and A. I. Larkin, Sov. Phys. JETP 58, 444 (1983). 

K. Frahm, Phys. Rev. Lett. 74, 4706 (1995); P. Brouwer and K. Frahm, Phys. Rev. B53, 
1490 (1996); B. Rejaei, Phys. Rev. B53, 13235 (1996). 

O. N. Dorokhov, JETP Lett. 36, 318 (1982); P. A. Mello, P. Pereyra and N. Kumar, 
Ann. Phys. (N.Y.) 181, 290 (1988). 

C. W. J. Beenakker and B. Rejaei, Phys. Rev. Lett. 71, 3689 (1993). 
V. Plerou and Z. Wang, Phys. Rev. B58, 1967 (1998). 

B. Jovanovic and Z. Wang, Phys. Rev. Lett. 81, 2771 (1998); K. Slevin and T. Ohtsuki, 
Phys. Rev. Lett. 78, 4083 (1997); C. Soukoulis et al, Phys. Rev. Lett. 82, 668 (1999); 
T. Ohtsuki et al, eprint |cond-mat/ 980922 1| . 

R. Landauer, IBM J. Res. Dev. 1, 223 (1957); D. S. Fisher and P. A. Lee, Phys. Rev. B 
23, 6851 (1981); E.N. Economou and CM. Soukoulis, Phys. Rev. Lett. 46, 618 (1981). 
K.A. Muttalib, J-L. Pichard, and A.D. Stone, Phys. Rev. Lett. 59, 2475 (1987). 
J-L. Pichard et al, J. Phys. France 51, 587 (1990). 

B. L. Altshuler, V. E. Kravtsov and I. V. Lerner, Sov. Phys. JETP 64, 1352 (1986) and 
Phys. Lett. A134, 488 (1989). 
[17] B. Shapiro, Phys. Rev. Lett. 65, 1510 (1990). 



7 



FIGURE CAPTIONS: 

Figure 1: Log-normal distribution of conductance P[ln(g)} given by eq. (1) in the in- 
sulating regime for two strengths of disorder, T = 0.2 (dashed line) and T = 0.7 (solid 
line). 

Figure 2: Distribution of conductance P(g) on the insulating and metallic sides of the 
crossover regime (r ~ 1) for two strengths of disorder, Y = 0.7 (solid line) and T = 1.6 
(dashed line). 
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